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c-j ■ Abstract 

q_j The soft function relevant for the production of an electroweak boson (7, W, Z or H) 

with large transverse momentum at a hadron collider is computed at next-to-next-to- 
leading order. This is the first two-loop computation of a soft function involving three 
light-cone directions. With the result, the threshold resummation for these processes 
can now be performed at next-to-next-to-next-to-leading logarithmic accuracy. 



1 Introduction 



The production of an electroweak boson, followed by its decay to leptons, is one of the most 
basic hard-scattering processes at hadron colliders. With more than 51b" 1 of data, the LHC 
experiments have now recorded millions of Z and W bosons decaying into lepton pairs, which 
allows for precision measurements, even at large transverse momentum of the boson. Given 
that Z's and W s at high px provide an important background to new physics searches and are 
used to calibrate jet energy scales, it is important to have good theoretical control of the cross 
section in this region. The complete 0(a 2 s ) corrections to vector boson production are known 
and have been implemented into numerical codes which allow for arbitrary cuts on the final 
state [H El [31 H] . However, since the p^-spectrum starts at 0(a s ), these codes only provide the 
next-to-leading-order (NLO) corrections for boson production at large transverse momentum. 
There is an ongoing effort to evaluate electroweak boson production in association with a jet 
to NNLO. The necessary two-loop results for the production of photons, Ws and Z's are 
known for some time [5J, El U\- Very recently, also the two- loop results for Higgs production 
in association with a jet were presented [5]. To obtain the transverse momentum spectrum 
at NNLO, these virtual corrections have to be combined with the real-emission corrections, 
which is a difficult task because of the very singular nature of the individual contributions. 

While the real-emission corrections are complicated in general, they simplify close to the 
partonic threshold. In this region, the hadronic final state consists of a single, low-mass jet 
and all radiation must either be soft, or collinear to the jet or the incoming hadrons. In this 
situation, the partonic cross section factorizes into a hard function times the convolution of a 
jet function with a soft function: 

a = H-J®S. (1) 

This factorization holds channel by channel, but with different hard, jet and soft functions. 
For 7, Z or W production, there are two channels, the annihilation channel q + q — > g + V 
and the Compton process q + g — > q + V. For Higgs production, the two relevant channels 
are g + g — >■ g + H and q + g — > q + H . A detailed derivation of the factorization formula in 
Soft- Collinear Effective Theory (SCET) [HI QUI ELT] and one-loop results for various ingredients 
were given in [12] . By now, almost all of the ingredients to obtain the threshold cross section 
to NNLO are known. The hard function H contains the virtual corrections, and can be 
obtained from the results in [5J El [8]. The two- loop jet functions J for quark and gluon jets 
were evaluated in [T51 [Ti] . The only missing ingredient is the NNLO soft function S which is 
computed in the present paper. 

The scale associated with the soft radiation is lower than the scale pt which is relevant 
for the hard function. As a consequence, the cross section contains perturbative logarithms of 
scale ratios which should be resummed. At next-to-leading logarithmic (NLL) accuracy, this 
was achieved in [T51 [121 ELY] . In the papers [121 EEB1 EES] the resummation was recently extended 
to NNLL. Using the results of [201 [2H [22], all the anomalous dimensions necessary for N 3 LL 
resummation were derived in [12]. With the soft function computed in the present paper, the 
resummation can thus now be performed at N 3 LL accuracy. 

The soft emissions are described by Wilson lines along the directions nf = Pi/Ei of the 
partons which participate in the hard-scattering process. To be able to derive results indepen- 
dent of the color representation of the partons, we use the color-space formalism introduced in 
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qg -> qV 
qg ->• qH 



gg -> gH 



Figure 1: The three Wilson-line configurations relevant for electroweak boson production. 

[23| [23] . The Wilson line for a particle in a color representation with generators T® is defined 
as the path-ordered exponential 

Si(x) = P exp (ig s j dtrii- A a {x + t m) lA . (2) 

If the particle is a gluon carrying a color index c, we have {T a )b c = —i fabc, and for an outgoing 
quark (or incoming anti-quark) with index f3 the generator is (T a ) a g = t^g. For an incoming 
quark (or outgoing anti-quark) with index /3, the generator is (T a ) a g = —t% a , which translates 
into anti-path ordering in (j2J) so that Sg(x) = S' q (x). 

The specific Wilson line configurations relevant for electroweak boson production in the 
different partonic channels are shown in Figure [TJ We will compute the soft function for the 
general case, where the initial-state Wilson lines are in representations 1 and 2, while the 
final-state Wilson line is in the representation J. We denote the directions of the associated 
Wilson lines by n\, n-i and rij, respectively. The soft function then has the form 

S(u) = J2(°\ S l S l &j\X){X\Si S 2 Sj\0) S(u - nj ■ px ) , (3) 

x 

with Si = Si(0). It measures the probability for soft emissions from the three Wilson lines 
with given momentum component u = nj ■ px along the jet direction. This component is 
relevant since the invariant mass of the final-state jet is given by the momentum p c ps Ej rij 
of the particles collinear to the jet and the momentum px of the soft particles as 

M 2 j = (p c + px) 2 ~p 2 c + 2Ejnj-p x . (4) 

The color indices of the Wilson lines Si are contracted with the color indices of the hard- 
scattering amplitude, while the conjugate Wilson lines Sj are contracted with the indices of 
the complex conjugate amplitude. These amplitudes collect the virtual corrections to the scat- 
tering process and their product is the hard function H in ([T]). However, the hard-scattering 
process involves only three partons and the color algebra of these amplitudes is diagonal [23] . 
As a consequence, the soft function reduces to a number times the trivial color structure 
S(u) = S(uj) 1 when applied to the hard function in a given channel. The function S(uj) 
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depends on the channel through the Casimir invariants associated with the representations of 
the three Wilson lines. 

By now, a number of NNLO results for soft functions involving two Wilson lines can be 
found in the literature, many of them appeared during the past year [251 127J [281 [291 USB EI] • 
For the case of multiple Wilson lines, only the anomalous dimensions were analyzed, first in 
the massless case [521 E3] and more recently also in the massive case [Ml EH EE]- In the 
massless case, the two-loop anomalous dimension was found to be proportional to the one- 
loop result. What came as a surprise at the time, is now understood to be a consequence of 
the strong constraints to which soft anomalous dimensions are subject [201 EH [221 EH l38| [39] . 
In our paper, we present the first full two- loop result for a soft function with three Wilson 
lines. Interestingly, these constraints, in particular the invariance of Wilson lines under a 
rescaling of the reference vector, drastically simplify our computation. As we will show below, 
all diagrams with attachments to the jet Wilson line vanish and the only non-zero diagrams 
are the ones which are present already in the case with two Wilson lines. 

In the next section, we turn to the evaluation of the diagrams. We first demonstrate that 
the soft function is color diagonal and that only a very limited set of diagrams contribute. 
We then discuss how the corresponding integrals can be evaluated and give the bare result for 
the soft function. In Section [31 we perform the renormalization and present our final result. 
Technical details and explicit results for all loop integrals needed in the computation are given 
in three appendices. 



2 NNLO calculation 

Before proceeding to the computation of the diagrams, we now demonstrate that the color 
structure of the soft function is diagonal and that the dependence on the light-cone vectors is 
completely determined by rescaling invariance. The color indices of the soft Wilson lines are 
contracted with the hard-scattering amplitudes. In color-space notation, the hard amplitude 
for the annihilation channel can be written in the form 

\Cqq-^gV )color = tp a Cqq^gY (5) 

where a, (3 and a are the color indices of the quark, anti-quark and gluon, respectively, and 
Cqq^gv is a function of the momenta and spins only. We have added a subscript to the vector 
\Cqq^ g v) to distinguish the color states of the hard function from the state vectors of the 
Hilbert space of the soft partons. The relevant color structure for Higgs production in the 
gluon channel is 

\Cgg-Jf.gH ) Color = ^fabc Cgg^,g}{ . (6) 

On general grounds, also the color structure d a bc could arise, but at least up to two-loop order 
this does not happen [8]. 

The diagrams contributing to the soft function to NNLO are shown in Figure [2J To 
understand why the color structure is diagonal, let us work out the color factor of the diagram 
D\ in the color-space formalism. The color structure associated with this gluon exchange is 
— Ti • T 2 = — T^T£. The minus sign arises, because one of the Wilson lines is conjugated. 
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Figure 2: Feynman diagrams that contribute to the soft function up to NNLO. In addition 
there are mirror symmetrical graphs, which we take into account by multiplying each diagram 
Di with a symmetry factor f h where fx = f 2 = f 3 = /s = 2, f 4 = 1 and f 6 = f 7 = 4. 
Additional diagrams, in which gluons attach to the jet Wilson line, vanish, see text. 

Its matrix element can be simplified as follows 

( 1' 2' J' | (-2\ • T 2 ) | 1 2 J) color = (1' 2' J' | 2\ • H + 2\ • Tj \ 1 2 J) color 

= d ( 1' 2' J'| 1 |1 2 J) color + ( 1' 2' J' \T X -Tj\12 J) color 

= ~ (Ci + C 2 - Cj) ( 1' 2' J' 1 1 11 2 J ) color , (7) 

where we have indicated with primes the different color state of the conjugate hard amplitude. 
In the first line, color conservation ^\ T" 1 = was used. The Casimir operator in the second 
line is T\ • T\ = C\\. The third line follows after applying color conservation two more times. 
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(a) (b) (c) 

Figure 3: Additional diagrams that turn out to be zero. Graphs (a) and (b) vanish because 
the associated integrals are scaleless, (c) because of its color structure. 

The color structure is thus trivial and it is convenient to define 

C s = 1 -(C 1 + C 2 - Cj ) = { Cf - Ca/2 for <^> (8) 
2 y C A /2 for qg -»■ q and gg -»■ g . 

Using color conservation and the commutation relations [T" 1 , Tf] = if abc Tf, the color factor 
of all diagrams can be expressed in terms of the Casimir invariants multiplying the unit matrix 
and the soft function takes the form 

( 1' 2' J'| S(u) |1 2 J) color = S(u)(l' 2' J'\ 1 11 2 J) color . (9) 

The color structure in the factorization theorem ([1]) is trivial and it can be written in terms 
of scalar hard, jet and soft functions. 

Not only the color structure, but also the dependence of the soft function on the light- 
cone vectors rij is very simple, because of the invariance of the Wilson lines under a rescaling 
rii — > A iii of the reference vectors. From the definition ([3]) it follows that under a simultaneous 
rescaling 

ni-)-Aini, 77,2-^2^2, nj->\jnj, u^-Xju (10) 

the soft function transforms as 

S(u)^±S(u). (11) 
The dependence of the soft function on the reference vectors must thus be of the form 



S(u) = — f(u) with Co = u \ — , (12) 

uj \ n x ■ nj n 2 ■ nj 

where the quantity Cj is invariant under the rescaling ( II Op . The factor of \/2 in Cj has been 
inserted for convenience. 

Having discussed its general structure, we now turn to the evaluation of the soft function to 
NNLO. Since the soft Lagrangian of SCET is identical to the ordinary QCD Lagrangian, the 
computation can be performed in QCD itself. At tree level only the vacuum state contributes, 
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and the soft function is trivially given by a delta-function. At NLO both virtual corrections 
to the vacuum state and one-gluon emission diagrams arise. In dimensional regularization 
the virtual graphs are scaleless and vanish. As the kinematic dependence of the soft function 
involves all of the reference vectors n,, see (I12p . the diagrams must have at least one gluon 
attachment to the Wilson line Si and one to S 2 - It follows that only the diagram D\ in Figure 
Ogives a non-vanishing contribution at NLO, while diagrams such as (a) and (b) in Figure 
|3] are scaleless and vanish. The calculation of diagram D\ reduces to the computation of the 
integral 

I x = I d d k 5(k 2 ) 6(k°) Ul ' n2 - 5(u -nj-k). (13) 
J rii ■ k n 2 ■ k 

This integral has been evaluated in [12] in an expansion in e = (4 — d)/2 up to the finite term, 
but here we also need higher order terms. As detailed in Appendix |A] we have computed this 
integral in closed form and obtain 

uj T(l — 2e) 

which is in agreement with the result from |12j . 

At NNLO the purely virtual corrections are again scaleless and vanish. Among the mixed 
virtual-real and the double real emissions, only the diagrams in Figure [2] are non-zero. The 
respective color factors can be worked out in the color-space formalism, as discussed above. 
We find C s for diagram D 2 , C 2 S for diagram D 3 , C S (C S — Ca/2) for diagrams D A and D 5 
and C s Ca/2 for the non-abelian diagrams D$ and D-j. Note that the diagrams in Figure |2] 
only involve attachments to the Wilson lines Si and S 2 . In addition there are diagrams with 
attachments to the Wilson line Sj. An example of such a diagram is the graph (c) in Figure 
El The color structure of this diagram is 

f abc T?T b 2 T c j, (15) 

which, by color conservation, vanishes when acting on three-parton states |1 2 J) co ior- In fact, 
the diagram vanishes in general: it could only be non-zero if it contained an anti-symmetric 
function of the three vectors rii, n 2 and nj. However, the only invariant Co is symmetric under 
the interchange ri\ n 2 and thus all diagrams with the color structure (Il5p vanish. The 
remaining abelian diagrams with attachments to all three Wilson lines also vanish, because at 
least one of the subdiagrams involves only two light-cone vectors and these integrals are thus 
scaleless. 

The technical aspects of the NNLO calculation are given in Appendix [Bj where we discuss 
the evaluation and give results for all integrals associated with the diagrams in Figure [2J After 
combining all contributions, we obtain the following result for the bare soft function 

s - H ^ H + ^I(|) 2 ' G , Sl(e) 



(Att) 2 cu \oj 



h ° s \° aSa ^ + V> 5 /( e ) + C sSs(e)] , (16) 
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where a s refers to the MS coupling constant, which is related to the bare coupling constant 
a° s via Z a a s fi 2e = e- e7B (4vr) e a° with Z a = 1 - /3 a s /(4vre) and (3 = 11/3 Cx- 4/3 7>n/. The 
variable u includes the dependence on the light-cone vectors and was defined in ( TT2l) . The 
one-loop coefficient reads 

Sl(f) ._ 8e ™.£W£!(L_f). (17) 

The three leading terms in e were given in jl2j. The two- loop coefficients are found to be 

_ , , 44 1 /4vr 2 268\ 22tt 2 1616 r „. 
^ (e)= -3? + 7(---J + ^-^7- + 56C3 

/ 9712 134tt 2 2728C 3 88tt 4 \ ml 2 , 



A 7 3e 2 9e 3 27 



/2624 _ 40tt 2 992C 3 \ 
V~81 9 9" J 



16 80 8tt 2 448 /2624 40vr 2 992C 3 \ ml 2^ 

-J e + C(e 2 ) 

s , (£) = _| + !|! + I^ + ^ £+cV ). (18) 



The function S haie (cj) is a distribution in u whose explicit form is obtained after expanding 



uj \ a I ne ml 



ln m (^ 



(19) 



The star-distributions are generalizations of plus-distributions to dimensionful variables, their 
definition can be found in [ID] , Because the expansion ffl9|) starts at 1/e, we evaluated the 
two-loop coefficients to 0(e). 

3 Renormalization 

The renormalization of the soft function is conveniently discussed in Laplace space. We intro- 
duce the Laplace transformed soft function 

f°° 1 
s(L,/i)= / due'"" S(u,fj,) with v= — — , (20) 

where, for later convenience, we have written s as a function of the logarithm of the Laplace- 
space variable 

L = \^ = \J-\\- ^'^ ) ■ (21) 
The Laplace transform of the bare function is immediately obtained using the relation 

i 

du e~™ or 1 "™ 6 = e~ ne7E r(-ne) n~ ne . (22) 
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After expanding in e, the Laplace transformed function is a polynomial in L. The function 1$ 
is also what is needed to perform soft-gluon resummation in the momentum space formalism 
of jH]. It fulfils the renormalization group (RG) equation 



dln/i 



s (L,fi) 



-AT 



cusp 



In --2T 



s (L,fi) 



(23) 



where r cusp denotes the cusp anomalous dimension, which is known at the three-loop level [12] . 
The anomalous dimension 7 5 was inferred to three loops in [12] using RG invariance of the di- 
rect photon-production cross section and the three-loop results for the hard anomalous dimen- 
sions [2"0"j I2"2"]. the quark jet anomalous dimension [13] and the Casimir scaling property of the 
soft anomalous dimension. Expanding the anomalous dimensions as r cusp = Yl^Lo ^n(f^) n+1 
and 7 5 = Yl^Loln (f^) n+1 5 one can easily solve the RG equation ( 1231) . To two-loop order, the 
solution takes the form 



s(L t /*) = 1 + (g) [2r L 2 + 2 7o s L + cf ] + 



An) 



2YlL A - 



ATr 



(A, - 3 7o s ) L 3 



+ 2 (l\ + ( 7o 5 ) 2 - /3o7o + r cf) L 2 + 2 ( 7 f + 7o 5 cf - /3 cf ) L + cf 
The expansion coefficients of the anomalous dimensions are 
r = AC S , Y x = AC, 

V y 

808 IItt 2 



Ca (? - t) ~ 2 i Tpn i 



(24) 



(25) 



7o 5 = 



H = C a C A \2K 



27 
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C s n f T F 



/224 4tt 2 



V 27 



9 



and the one- loop constant is cf = C s ii 2 . The color factor C s depends on the partonic channel 
and was defined in (jH]). 

The function s renormalizes multiplicatively, s = Z s s bare ; and Z s fulfils the same RG 
equation (1231) as the renormalized soft function. Solving this equation [201 [22], one derives the 
following expression for the logarithm of the Z-factor, 



In Z Q 



An 



[ An 



(2r L + 7o 5 ) 



— (2r L + 7o ) - — + _(2r 1 L + Tl ) 



4e 3 2e 2 



(26) 



Since all the necessary anomalous dimensions are known, the Z-factor is completely deter- 
mined at the two-loop level. The cancellation of all divergences l/e n , for n = 1...4, in 
the renormalized result provides a strong check of our calculation. We finally obtain for the 
non-logarithmic two-loop coefficient 



_s_ r2 vr 4 Z2428 3357T 2 22C 3 14tt 4 \ / 656 50vr 2 8( 3 \ 
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which is the main result of our paper. Note that the coefficient of the color structure C 2 is 
just one half of the one- loop coefficient squared. This is a consequence of the non-abelian 
exponentiation theorem for Wilson lines. Note that the relation would be more complicated 
for the momentum space function. Since the position-space and Laplace-space functions are 
related by Wick rotation, the non-abelian exponentiation takes the same form as in position 
space. 

To illustrate the size of the corrections, we now evaluate the soft function numerically for 
rif = 5. For qq — > gV, where C s = Cf — Ca/2, we find 

s(L,n) = 1 + (Jp) (-1.333L 2 - 1.645) 

/ QJ \ 2 

+ \T) ( - 889l4 + 6 - 815l3 - 7 - 018l2 + 6 - 17L + 13 - 22 ) > ( 28 ) 

while the case C s = Ca/2, which is relevant for qg — > qV and for Higgs production, yields 

s(L,//) = l+(^) (12L 2 + 14.80) 

+ {^j 2 (72L 4 - 61.33L 3 + 260.6L 2 - 55.53L + 2.757) . (29) 

In the first case, the corrections are quite small: setting a s = 0.1 and varying —2 < L < 2 
the two-loop contribution is below half a per-cent. The situation is somewhat different for 
C s = Ca/2- While the constant piece is small, the corrections can be sizable, if the soft 
function is not evaluated at its natural scale: for a s = 0.1 and L = —2, for example, the 
perturbative expansion reads s(L, fj) = 1 + 0.50 + 0.18. In a resummed computation, the 
hard, jet and soft functions are evaluated near their natural scale and so we do not expect 
large corrections to the results obtained in [18]. However, to make definite statements, the 
different two-loop ingredients should be combined, which will be done elsewhere. 

4 Conclusions 

We have computed the soft emissions for electroweak boson production at large transverse 
momentum near the partonic threshold to two-loop order. The relevant soft function involves 
light-like Wilson lines in three directions: along the beam directions and along the direction of 
the final-state jet recoiling against the electroweak boson. The invariance of the Wilson lines 
under a rescaling of the reference vectors puts strong constraints on the form of the result 
and greatly simplifies the computation. In particular, the diagrams which involve emissions 
or absorptions from the jet Wilson line vanish in dimensional regularization, such that the 
computation ends up being similar to earlier computations of dijet soft functions. 

With our result, the last remaining ingredient to extend the threshold resummation for 
electroweak boson production processes to N 3 LL is now available. In the future, we will 
combine our result with the two-loop jet and hard functions to obtain improved predictions 
for transverse momentum spectra. 
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A Evaluation of the NLO integral 

In this appendix we outline the calculation of the NLO integral 

I x = f d d k 6(k 2 ) 9(k°) n \' n2 , 5(oj -nj-k). (30) 

J Hi ■ k 77-2 • k 

We closely follow the strategy adopted in [T2]. First we introduce light-cone coordinates and 
decompose the four-vectors as 

\/2n x ■ n 2 y/2ni ■ n 2 

with n 2 = and rij • q± = 0. The result for the integral is frame independent, but it is 
convenient to work in a reference frame where n\ and n 2 are back-to-back. Notice that this 
does not imply n\ ■ n 2 = 2. We may, for instance, choose a frame where ri\ = (1,0, 0,1) and 
n 2 = (A, 0, 0, —A) with A = n\ ■ n 2 /2. In light-cone coordinates the integral takes the form 

//*oo rl 
dk + dk_ J d\k ± \\k ± \ d - 3 J dcos9sm d - 5 9 

x - kl) 0(kP) ~ + y/rwjl$±\ cos 6>) , (32) 

where 9 is the angle between k± and nj± and Q n is the n-dimensional solid angle. We next 
use the on-shell condition to perform the [-integration. Parameterizing 

k + = — ^- xy, k- = — ^- (1 — x)y , (33) 

nj- n J+ 

we arrive at 

r\ pi poo pi 

h = u~ 2e dx dy (xx)- 1 -" y' 1 ' 2 " / d cos 9 sin d " 5 9 6(1 - y + 2y^fx~x~ cos 9) 
2w Jq J J_! 

(34) 



with x = 1 — x and Q as defined in ffT2J) . We then use the delta-constraint to perform the 
^-integration. The subsequent integration over cos 9 gives rise to an hypergeometric function 
and one is left with 

h = fffTej I I dx ( x5; )" 1 " e 2Fl (I ~ e ' ~ e 5 1 " e ' Ax ^) ■ ( 35 ) 

After performing a variable transformation to the argument of the hypergeometric function, 
the integral can be solved in closed form and we obtain 

h = -— U T(l-2e) ■ (36) 
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B Details of the NNLO calculation 



The calculation of the NNLO diagrams proceeds along the same lines as the NLO calculation 
that we outlined in Appendix |A] Here we briefly describe the new elements that arise at 
NNLO, and give results for the individual integrals associated with the diagrams in Figure [2J 
We start with diagram D 2 . As the soft function only depends on the total momentum of the 
radiation, it is convenient to express this diagram through the cut of the vacuum polarization 

U^(k) = (k 2 9ilv - kpK) U(k 2 ) . (37) 

Diagram D 2 then involves an integral of the form 



The contribution from k^k v vanishes, since the light-cone propagators cancel out in this term 
and one is left with a scaleless integral. Inserting the one-loop expression for the vacuum 
polarization, we obtain the integral 

u = / a- »<*•) ni . k 27 ik ^, *<« - ■ *) ■ (39) 

The calculation of this integral is similar to the one sketched in Appendix IAJ except that we 
now use a slightly different parameterization given that k 2 ^ 0, 



k + — xy, k_ = (1 — x)y , \k±\ — \/(l—x)xyu. (40) 

n.j- n J+ y/nj+nj- 

We then use the delta-constraint to perform the ^-integration. The remaining integrals are 
standard and can be performed along the same lines as in Appendix |A1 yielding 

to r(l-4e) V ; 

Next, we turn to the evaluation of the abelian diagrams -D3, and -D5. The planar 
diagram D3 involves the integral 

[<*] / w , ( ;"'", 2>2 i| ";° r , t ,"T" „ ■ < 42 > 

L J J L J m • k (m ■ k + m ■ I) n 2 ■ I (n 2 ■ k + n 2 ■ I) 

where we introduced the short-hand notation [dk] = d d k 5(k 2 ) 9(k°). Using partial fractioning 
identities this integral can be expressed as I 3 = J 4 /2 — I 5 , where 7 4 and I 5 are integrals 
associated with diagrams -D4 and D§. Their precise definition will be given below. The 
integral in diagram D4, 

[dk] [ [dl] {Ul ■ n2) ] 6{U ~ nj \ k ~ W / • , (43) 
J rii ■ k rii ■ I n 2 ■ k n 2 ■ I 
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can easily be calculated since it is just the convolution of two NLO integrals, 

' du' h(uj - <J) h(J) = r2(1 + £) ,7~ £) . (44) 

The non-planar diagram D 5 is more complicated. It gives rise to the integral 

[dk] [ [di\ ^ i - n f 6iUJ - nj \ k ( - nj : l) n . (45) 

As the delta-function only constraints the sum of the cut momenta, we first combine them to 
q = k + 1 and rewrite the integral as 

, d {nx • n 2 ) S(lo - nj • q) f f (m ■ n 2 ) S d (q - k - I) 



d d q V1 z/ v J ^ / [dk] / [dl] v 1 z/ - 7 ^ ' . (46) 

rif qn 2 - q J J n x -kn 2 -k 

The integrations over k and I can be performed using the auxiliary integral 058)) from Appendix 
ICl We further adopt the parameterization (1401) together with an integral representation of the 
hypergeometric function and arrive at 



n d _ 3 . rri-e) 



dx dy {xx)~ l - 2e y- 1 - 4 * / dcos9sm d - 5 e 
Jo Jo J-l 



,l+4e l^zj V 



lo T(l - 2e) 7 

y cfo m 1 " 26 (1 - it 2 )" 1 - 26 v- 1 - 6 (1 - ^ 2 ) e <?(l-| + j/Va^ucos0) . (47) 

We now use the delta-constraint to perform the y-integration. The subsequent integrations 
over cos 8, x and v are standard and can be performed along the same lines as in Appendix IA1 
This yields an integral over a product of two hypergeometric functions, which we solve order 
by order in the e-expansion. The result reads 

^ ~-4e T 2 (l-6) (I QA 27T* , 



^ " ' -T-8C 3 -^ + 0(e 2 ) • (4 



u T 2 (l-2e) \e 

We finally turn to the non-abelian diagrams D e and D 7 . For the integral associated with 
diagram D 6 , 

J J rit- k (m- fc + 711-Z) K- k + n 2 -l) (k + I) 2 ' v ; 

we again combine the cut momenta to q — k + I and make use of the auxiliary integrals 
from Appendix O The subsequent steps are by now straightforward and we obtain a closed 
expression 



7T 



2-2<e 



..„r(i + 2 t) r(-2^H^ 



u (1 - 2e)r(l - 4e) 
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The last diagram D 7 contains both a one-particle and a two-particle cut. The one-particle cut 
involves the real part of the loop integral 

d d k 



(ni • q + n\ ■ k) n 2 ■ k (q + k) 2 k 2 

2- £ ^r 2 (i + e)r 3 (-e) ( m-n 2 



= 2ir 2 - e e~ zne v — ' , v — '- , (51) 

r(-2e) \2n x -qn 2 -q) ' 1 ; 

which gives rise to 

hi = [ d d q 5(q 2 ) 9(q°) ( Wl ' " 2 ) ^ 5(u - nj ■ q) . (52) 
J \2ni-qn 2 -qJ 

This integral is similar to the NLO integral that we discussed in detail in Appendix |A1 It can 
be computed along the same lines and gives 

hl - ^7 U T(-Ae)T 2 (l + e) " (53) 
Finally, the two-particle cut of diagram D7 leads to the integral 

t fun f un (ni-n 2 ) (n 2 - k + 2n 2 - 1) 5(u -nj ■ k -nj -l) 

For this integral we follow the same strategy that we adopted for the integral I5. We again 
encounter a product of two hypergeometric functions, which we solve order by order in the 
e-expansion. The final result reads 

2vr 2 - 2e . T 2 (l - e) / 1 7T 2 41 5tt 4 ,A , x 

h.2 = ui~ Ae n \ -r - — - -Ca - -5-6 + 0(e 2 ) . (55) 



u T 2 (l-2e) \e 3 4e 2 

C Auxiliary integrals 

For the evaluation of the diagrams with two-particle cuts, the following integrals are useful: 
[dk] J [dl] 5 d (q -k-l) = vr 1 - 2 ^ ~ ^ (g 2 )" e , (56) 

w/H '-,t5,-.r o 

.l- e r (~ e ) f 2n l ■ V ™2 • qY f 2 \-l-2e r^f .. K " n a) |<f± | 2 



r(l - 2e) \ ni • n 2 J \ 2n 1 -qn 2 -qJ 
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J [dk] J [dl] 



(m ■ n 2 ) 5 d (q — k — I) 



= 7T 



1-e 



2 Til • k 7l2 ■ I 

r(-e) (2n x -qn 2 -q 



) 



{q 2 )- 6 2F 1 -e 



e; 1 - e; 



(ni ■ n 2 ) q 2 




r(l-2e) V wi - "a 



2 ni • q n 2 • q 



where [dk] = d d k 5(k 2 ) 8(k°). Note that the above relations imply q 2 > and q° > 0, since 
the momentum q is the sum of two physical momenta. Our results are in agreement with [30] . 
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